Vortex solitons in superconducting films with periodically modulated thickness 
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We report magnetoresistance measurements through Nb films with a periodic thickness modula- 
tion. The pillar shaped large-thickness regions of the sample, which form a square lattice, act as 
repulsive centers for the superconducting vortices. For low driving currents along one of the axes 
of the square lattice, the resistance R increases monotonously with increasing magnetic field B and 
the R-B characteristics are approximately piecewise linear. The linear R vs B segments change 
their slope at the matching fields where the number of vortices is an integer times the number of 
pillars in the sample. Numerical simulations allow us to associate the different segments of linear 
magnetoresistance to different vortex-flow regimes. We argue that the magnetoresistance below 
the first matching field is governed by single-vortex physics while at higher fields is dominated by 
different vortex-soliton regimes. 



I. INTRODUCTION 

The dynamics of driven superconducting vortices has 
been the subject of intense theoretical and experimental 
research during the last decade. This has been mainly 
fueled by the technological interest in reducing the dissi- 
pation associated to the vortex motion and by the pos- 
sibility of using vortex matter as a model system with 
tunable parameters to study static and dynamic phases 
and their transitions. In vortex systems it is possible 
to control the force driving the vortices by applying an 
electric current across the sample, to set the vortex den- 
sity through external magnetic fields, and to tailor the 
potential landscape using lithographic techniques. 

The fact that the density of vortices in a sample can be 
tuned applying an external magnetic field, has allowed a 
number of studies that analyze commensurability effects 
between the vortex lattice and an underlying periodic 
potential. The typical matching signature is the appear- 
ance of minima in the resistance (or maxima in the crit- 
ical current) for magnetic fields such that the number of 
vortices riy is a multiple (or a simple fraction like 1/4, 
1/2 and 3/4) of the number of minima rip in potential 
landscape. For this purpose, films having different typ es 
of artificial pinning centers like holes P™» magnetic^ or 
non-magnetic particles^ in different geometries^^ have 
been studied both experimentally and theoreticallyPHHI 

Driven vortices in periodic potentials present a rich va- 
riety of phenomena, 28 including transitions between dif- 
ferent vortex flow regimes. A prominent example is the 
rich phase diagram predicted by Reichhardt et alW^ for a 
square pinning potential by means of numerical simula- 
tions, which was later confirmed experimentally (see Ref. 

A common feature in the general problem of inter- 
acting particles in a periodic potential is the formation 
of incommensurate structures! 3 ^ These appear as a re- 
sult of the competition between inter-particle interac- 
tions, which favor a specific crystal structure, and the 



external potential, which in general, favors another. In 
vortex systems these discomme nsura tions have been pre- 
dicted in numerical simulation d 31 ! 32 ! and also visualized 
experimentally,^ but their dynamics and the experimen- 
tal consequences of their presence have been less ex- 
plored. 

In this work we analyze the dynamics of driven vor- 
tices in a square pinning potential as a function of the 
vortex density, the driving force, and the lattice param- 
eter of the pinning potential. To that aim, we generate 
a periodic thickness modulation in Nb films and perform 
magnetoresistance measurements in the superconducting 
state. The potential energy of a vortex increases with 
the sample thickness due to the extra energy associated 
to the vortex lengthP32l! This allows to set the geome- 
try and the intensity of the pinning potential with high 
precision using lithographic techniques. 

For a film with a square lattice of protruding cylinders 
we observe clear commensurability signatures. The re- 
sistivity increases linearly with increasing magnetic field, 
but instead of the usual minima at the matching fields 
we observe changes in the slope of the R vs B curve. 
To interpret the results we perform molecular dynamics 
simulations, and argue that the observed behavior can be 
ascribed to weak vortex- vortex interactions (compared to 
the interaction with the periodic potential) and the pres- 
ence of moving discommensurations (solitons). 

The rest of the paper is organized as follows: In Sec. 
|Il| we describe the experimental setup. In Sec. |III| we 
present magnetoresistance measurements illustrating the 
commensuration effects. In Sec. |IV| we present the model 
and methods to describe the driven vortex dynamics. In 
Sec. [V] we present the results of numerical simulations 
that reproduce qualitatively the experimental observa- 
tions. In Sec. |VI| we present a simplified one-dimensional 
model that reproduces qualitatively the main features of 
the experimental results. Finally, in Sec. | VII| we present 
our concluding remarks. 
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FIG. 1. a) Scanning electron microscopy image of the mea- 
surement circuit. The current path (coloured in red) is cou- 
pled to contacts for voltage measurement (coloured in yellow) . 
The box of dotted style lines in the current path indicates the 
region where thickness modulation is performed, b) Scanning 
electron microscopy image of the square lattice of resist cylin- 
ders before performing the RIE. c) Image of a sample with 
Nb cylinders of diameter d — 350nm and lattice parameter 
a = 570nm. 



II. EXPERIMENT 

We fabricated a circuit (see Fig. [IJ consisting of an 
electrical current path 60um wide coupled to electrodes 
for longitudinal voltage measurement!-^ This structure 
is fabricated in two stages. First, the desired pattern for 
the extension electrodes is defined with optical lithogra- 
phy on a AZ-9260 resist layer. Nb is sputtered on top 
and a lift-off process is performed. Second, using the 
Nb deposited as an alignment mark, the pattern for the 
current path and the first part of the contacts is written 
with e-beam lithography on a bilayer of MMA / PMMA 
resists. Then again, a Nb film is sputtered on top and a 
lift-off process is performed. 

The periodic modulation of the thickness is generated 
on 60^m long regions located on the current path. Us- 
ing the negative photo-resist ma-N 2403 a square lattice 
of resist cylinders (see Fig. I]) is fabricated by e-beam 
lithography. When the sample is exposed to reactive ion 
etching (RIE), the cylinders protect the covered portion 
of the film, generating the desired thickness modulation. 
During the RIE process we work with a mixture of SFq 
and At at a pressure of 1.33Pa, and an acceleration volt- 
age for the ions of approximately 110U. With these pa- 
rameters we obtain a Nb etching speed of 80nm/min. 
After the RIE process, the sample is immersed four hours 



FIG. 2. Resistivity vs magnetic field in a Nb film with a square 
lattice of protruding cylinders at T = 7.52K (T/T c = 0.996) 
(sample I). The lattice parameter is a = 570nm {B^ = 62 G) 
and the diameter and height of the cylinders is d = 350nm 
and h cy i — 40nm, respectively. The resistance presents an 
approximately piecewise linear behavior with slope changes 
at the first two matching fields. As the current driving the 
vortices is increased (see inset) the kink at the second match- 
ing field B — 2£>0 becomes less pronounced and eventually 
disappears. 



in acetone in order to remove the resist. Figure [T] shows 
a scanning electron microscopy image of one sample with 
cylinders of diameter d = 350nm and lattice parameter 
a = 570nm. As usual, a determines the first matching 
field = 4*o /a 2 , where 4>o is the flux quantum. At 
B = B,f, there are as many vortices as unit cells in the 
cylinder lattice (i.e. ny =np). 

In what follows we will present magnetoresistance mea- 
surements for two samples: sample I with B$ — 62 G, 
and sample II in which B^ = 48 G. The film thickness is 
150 nm and 110 nm for samples I and II respectively. In 
both samples the nominal cylinder's height with respect 
to the film is h cy i = 40 nm. 



III. MAGNETORESISTANCE 
MEASUREMENTS 

The magnetoresistance measurements were carried out 
in a four probe geometry within a pulse tube refrigerator 
with 2.3K base temperature. The voltage was measured 
at constant temperature T and current J along the cur- 
rent path (see Fig. [T]) using a nanovoltmeter. The mag- 
netic field was increased from zero up to a maximum 
value of B max ~ 3-B^ and then decreased to zero again. 
No hysteresis effects were observed. 

In Fig. [2] and Fig. [3] we present magnetoresistance 
measurements for sample I. In what follows we focus our 
analysis on the data obtained within the range of temper- 
atures 7.5K < T < 7.54K < T c , where T c ~ 7.55K is the 
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FIG. 3. Resistivity vs magnetic field in a Nb film with a 
square lattice of protruding cylinders for J = 26.6 A /cm 2 and 
different values of the temperature (sample I). 



superconducting critical temperature of the sample, and 
the range of current densities 6 Aj cm 2 < J < 27 A/ cm 2 . 
For lower temperatures or currents the commcnsurabil- 
ity effects are weak, presumably due to the dominance 
of the intrinsic pinning mechanisms over the periodic 
potentialpSEH For temperatures or currents above the 
selected intervals, we observed a finite dissipation at zero 
external magnetic field which indicates the presence of 
vortices induced by the external current or other dissipa- 
tion mechanisms. 

The most important result is the presence of commen- 
surability effects at the first two matching fields B = B^, 
and B — 2B$ for the lowest temperatures and currents 
considered (see Fig. |2| . Instead of the usual minima we 
observe changes in the slope of the R vs B curve. This 
behavior can be attributed to vortex-vortex interactions 
whose intensity at distances of the order of the lattice 
spacing of the periodic potential is sufficiently weak com- 
pared to the interaction with the pinning potential. In 
order to observe the usual minima in the magnetoresis- 
tance, the vortex lattice must be rigid enough such that 
away from the matching fields the total force on the lat- 
tice due to the pinning potential averages to a small value 
or zero. In a rigid lattice, which is expected for strong 
vortex-vortex interactions, the pinning force on the lat- 
tice is determined averaging the pinning force over all 
vortices. Away form the matching fields, the force due 
to a periodic pinning potential has different directions 
for different vortices which leads to cancellations and a 
relatively small pinning force on the lattice. Near the 
matching fields, however, the force on many vortices of 
the lattice has the same direction and the contributions 
add up to increase the pinning force, producing a slow- 
down of the vortices and a minimum in the resistance. 

For low enough magnetic fields B ~ 0, the average 
inter-vortex distance (r) is large compared to the pene- 
tration depth A, and the typical vortex- vortex interaction 



is exponentially small cx '/\^In this regime, we ex- 
pect the resistance to be given by the number of vortices 
multiplied by the average velocity (v) sv of a single vor- 
tex when it travels across the sample. If (v) sv ^= wc 
expect a linear increase in the resistance with increasing 
magnetic field: 



R oc n 



v {v) 



(1) 



where ny = AB/<f>o is the number of vortices in the 
sample, A is the area of the film. 

In our measurements (see Fig. [2| we observe an ap- 
proximately linear behavior up to the first matching field 
i?0 where the number of vortices is equal to the number 
of minima of the pinning potential. This linear behavior 
and the absence of a minimum in the magnetoresistance 
at the first matching field indicate a weak effect of the 
vortex- vortex interactions for B < B^. For fields B > B<p 
we observe again a linear increase of the resistance with 
increasing B but with a different slope. The behavior of 
the resistance for B^ < B < 2B^ is well described by 

R(B cj} < B < 2B< i> ) oc up (v) sv + (n v - n P )v so i (2) 

where v so i > (v) sv . This result can be interpreted as 
follows: the additional vortices over the first matching 
field n so i = ny — np have an average velocity v so i , while 
the average velocity of the remaining np vortices is given 
by the single vortex value (v) vs . As we will see in Sec. 
|V| through numerical simulations and in Sec. V 
a toy model, the additional n so i vortices add n so i 



using 
vor- 
tex solitons to the system which have an average velocity 
v so i. Each additional vortex generates a discommensu- 
ration and produces a cascade of vortex displacements 
that moves faster than the individual vortices. There are 
no vortices overtaking others and all vortices that par- 
ticipate in the soliton at a given time have their average 
velocity increased. However, each soliton can be effec- 
tively viewed as a vortex moving with a larger velocity 
than the background vortices which maintain their aver- 
age velocity. 

The linear behavior of the resistance observed for 
B > B^ suggests a weak soliton-soliton interaction as 
the soliton velocity v so i is independent of the number of 
solitons. Interestingly, for the lowest temperatures and 
currents considered, in sample I we observe a linear in- 
crease in the resistance for fields B > B^ up to B = 2B^, 
where a second increase of the resistance of the form 

R(B > 2B<p) tx n P (v) sv + n P v sol + (n v - 2np)v 2so i (3) 



is observed. Again, a simple model (see Sec. VI) suggests 
that for B > 2B^ the extra vortices can generate mov- 
ing discommensurations (solitons) with average velocity 
V2soi > v soi i over a background of np solitons moving at 
velocity v so i. 

Upon increasing the external current (i.e. the vor- 
tex driving force) or the temperature (which increases 
the range of the vortex- vortex interactiona-E^l) the kink at 
2B& becomes less pronounced and eventually disappears 
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FIG. 4. Resistivity vs magnetic field in a Nb film with a 
square lattice of protruding cylinders at T — 7.65K (sample 
II). The lattice parameter is a — 650nm = 48 G) and the 
diameter and height of the cylinders is d = 340nm and h cy i = 
40nm, respectively. The resistance presents an approximately 
piecewise linear behavior with a slope changes only at the 
first matching field. As the current driving the vortices is 
increased, the slope of the low field segment increases, while 
it remains approximately constant for B > B^,. 



(see Figs. [3] and [4] ). In the parameter regime where 
the kink at IB^ is not observed, the slope of R{B) for 
B < B^ increases with J and T but it remains approxi- 
mately constant for B > B^. This implies that the single 
vortex average velocity (v) sv increases much faster than 
the soliton velocity v so i with increasing J and T. As 
we will show in the following sections, (v) sv is governed 
by the difference between the Lorentz force Fc and the 
maximum force generated by the periodic potential Fp, 
while in this regime v so i is to a good approximation de- 
termined by Fq alone. Under these conditions, a small 
relative change in the external force (i.e. on J ) can pro- 
duce a large relative change in (v) sv . An increase in the 
temperature produces a similar effect reducing the max- 
imum pinning force which leads to an increase in (v) sv 
while leaving v so i approximately unchanged. 



IV. MODEL AND METHODS 

Molecular dynamics has proven to be a very success- 
ful technique to describe the static and dynamic prop- 
erties of vorte x systems in a variety of experimental 
situations! 29 * 32 ^ Here we model the magnetoresistance ex- 
periments described in the previous section. We con- 
sider an effective two dimensional model to describe the 
vortices in the superconducting film. The overdamped 
equation of motion for the i-th vortex is given by 



where = r*i = (xi,yi) is the velocity, rj is the Bardeen 
Stephen friction, 



rvv \ ^ f 



(5) 



is the force due to the interaction with the other vortices, 
V(fi) is the external potential, and Fc = —Jxz, a driv- 
ing force generated by the transport density supercurrent 
J. 

In the limit A>(, we have 



fij = foKi 



13 



(6) 



where r,- 



■ y — ?i — rj, K\[r/X) is the modified Bessel func- 

tion of the first kind, and /o = g^rp- ■ Since Ki(r/\) oc 
exp(— r/A) for r > A, we cut off the vortex- vortex in- 
teraction for distances larger than 7. 5 A. In what follows 

we take A, — ^ and fo as distance, time and force units, 
Jo 

respectively. 

We model the interaction between the i-th vortex and 
the square lattice of cylinders with the potential: 



F P 



2vr 



. 2lTXi , 2-KVi 

sm( ) + sm( J 



(7) 



where a is the lattice parameter of the cylinders struc- 
ture, and Fp is a function of a and A.^ 

To simulate the magnetoresistance measurement we 
consider a system with periodic boundary conditions and 
an external force along one of the principal axes of the 
square lattice (the y direction). We start with Uq vor- 
tices, add at random positions An vortices and perform 
8000 molecular dynamics (MD) steps (At = 0.007) to 
reach a stationary state. We average the quantities of 
interest over 30000 additional MD steps. In particular 
we calculate the total average velocity (V) = *}™ v (vi) 
whose projection (V y ) along the y-axis is proportional to 
the measured voltage drop across the sample associated 
to the motion of the vortices. 

In the low density limit the behavior of (V y ) can be 
understood assuming that the inter- vortex interaction is 
negligible. We may then calculate the average velocity 
as (V v )ni — riv(v) sv where (v) sv is the mean velocity 
for an isolated vortex in the system, and determines the 
slope of (Vy). The motion of a single vortex in the poten- 
tial generated by the cylinders follows a one-dimensional 
path parallel to the y axis equidistant to two rows of 
cylinders. Along these paths, the pinning potential has 
the sinusoidal form 
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aF P 
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sin(27r?//a). 
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with an associated pinning force Fip>(y) = —dVip>(y) /dy. 
For Fc larger than the maximum pinning force Fp = 



5 



max[Fm(y)], the average velocity can be calculated in- 
tegrating the equation of motion: 



(v)sv = 



dy 

JO [l-F 1D (y)/F c \ 



Fc, 



(9) 



which can be evaluated for the potential of Eq. Q to 
obtain 



F 2 c 



F 2 



(10) 



Equation ([8| shows that (v) sv vanishes as y/Fc — Fp for 
Fc — Fp — > 0, which is a generic behavior for potentials 
having a quadratic maximum in the pinning force. 

To select the parameters in the numerical simulations 
we consider that: i) The linear behavior of the magne- 
toresistance at fields below the first matching field and 
the absence of a minimum at it indicate that Fc > Fp 
and weak vortex- vortex interactions for B < B^. ii) The 
large relative increase in dR/dB at the first matching 
field indicates that (v) sv is much slower than the maxi- 
mum velocity expected in the absence of a periodic po- 
tential v max = Fc/t], which implies that Fc — Fp <C Fe- 
in) The value of the penetration length for Nb at low 
temperatures^ is A(T — > 0) ^ 45 nm which is much 
smaller that the nearest neighbour cylinder spacing a in 
our films. It increases with increasing the temperature 
and diverges at the superconducting critical temperature 
T c . 

We will therefore focus from hereon on the weak inter- 
action A < a and weak driving force Fc — Fp regimes. 
The range of the interaction cannot be taken too short, 
however, because the single vortex behavior would extend 
to fields above B$ which is not what is experimentally 
observed. 



V. TWO-DIMENSIONAL SIMULATIONS 

In Fig. [5] we present numerical results for the total av- 
erage velocity along the y axis (V y ) vs ny/np = B/B<p for 
different values of the external force Fc- For forces over 
a threshold Fp — 2 and low vortex densities (ny < rip), 
the average velocity increases linearly with the number of 
vortices in the system following the expected behavior for 
non- interacting vortices (thin lines in Fig. [5| . The linear 
behavior of (V y ) persists up to the first matching field 
ny — np where there is a kink in the curves. In qualita- 
tive agreement with the experimental observations (see 
Fig. [4]), at low vortex densities there is a large relative in- 
crease (> 100%) in the slope of the curves for a moderate 
(< 5%) increase in Fc, while the slope remains relatively 
constant for ny > np. The large relative increase in 
(v) sv with increasing Fc can be easily understood from 
Eq. ([lO]) in the regime A F — F c - F P < F C - 

For ny > np, (V y ) presents again an approximately 
linear behavior. The deviation of (V y ) from the single 



2.5 



& 1.5 



1 - 



0.5 



9! 



3-e °-e-e-e>-3 &>-e-&4 Q-q-^o 



bbbB fa - a - aBBaaaE- ^ □ 








'Fc 
Fc 
F C 



2.02 
2.05 
2.09 



0.5 



ny/np 



1.5 



2.5 



vortex result for ny > np indicates a more active role of 



FIG. 5. Total average vortex velocity vs the number of vor- 
tices in a 15 x 15 square array of cylinders with lattice pa- 
rameter a = 6. For Lorentz forces Fc larger than Fp and 
ny/np < 1 the average velocity is well described assuming 
non-interacting vortices. At nv = np an increase in the slope 
is observed which is related to the propagation of solitons for 
ny > np. The dot-dashed style lines present the results ex- 
pected for non-interacting vortices (see text). For Fc = 2.02 
the inset contains velocity profiles along a line in the y direc- 
tion when there are 2 (□), 3 (o) and 4 (o) extra vortices on it. 
For clarity the curves are shifted 1.5 units vertically. In this 
regime, each additional vortex generates a disconmensuration 
that travels faster than the average vortex velocity. 



the vortex- vortex interactions in this regime. The inset 
in Fig. [5] shows velocity profiles along a line in the di- 
rection of Fc ■ It can be seen that each additional vortex 
on that line adds a disconmensuration that travels faster 
than the average velocity. As we argue in Section |LTl} the 
slope of R(B) in this regime is related with the propaga- 
tion of these structures. To understand this wc study the 
vortex velocity distribution P(V) for different densities. 
The single vortex probability distribution P SV (V) can be 
calculated exactly for the potential of Eq. ^ (see Sec. 
VI and Appendix |A|) . P SV (V) is peaked at Fc — Fp and 
Fc + Fp, the minimal and maximal velocities, respec- 
tively, for a single vortex in the periodic potential. The 
peak at Fc — Fp is more pronounced and for A^? <C Fp 
it dominates the average velocity. There is a large prob- 
ability of finding a vortex with a slow velocity because 
it spends more time 'climbing' the pinning potential bar- 
rier than 'descending' it. In Fig. [6] we plot P(V) for 
ny/np = 0.64 and 1.06. In the first case the numerical 
result is in very good agreement with whats expected for 
non interacting vortices, while in the second case, two 
new peaks appear. The inset in Fig. [6] shows the devi- 
ation of P(V) from the single vortex result. The proba- 
bility of finding a vortex with a slow velocity ~ Fc — Fp 
decreases, and this results in an increase on the average 
velocity. The two extra peaks in P(V) are associated 
to the velocities of a vortex pushing or being pushed by 
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that vortex-vortex interactions are negligible at large dis- 
tances. 

The overall behavior of the magnetoresistance depends 
on the intensity of the vortex-vortex interactions at dis- 
tances of the order of the lattice parameter a of the pe- 
riodic potential. For weak interactions |/ij(a)| <C Fp we 
expect the measured resistance at low densities (ny < 
np) to be given by Eq. 0, with a single vortex average 
velocity 



FIG. 6. Vortex velocity distribution P(V) vs (V - F C )/F P . 
For densities ny/np < 1, P{V) = P SV (V), while for ny jn P > 
1 there are two additional peaks (see arrows) and the prob- 
ability of finding a vortex with the slow velocity Fc — Fp 
decreases (see inset). 



another vortex in the moving discommensuration. 

The velocity v so i that can be extracted from the curves 
of Fig. [5] is ~ 80% of F c indicating that the discommcn- 
surations move with a velocity close to the one expected 
for the vortices in the absence of pinning potential. 



VI. ONE-DIMENSIONAL (ID) TOY MODEL 

In this section we consider a ID toy model that pro- 
vides a simple qualitative picture of the vortex dynam- 
ics in the experimentally observed regimes. The two- 
dimensional numerical simulations of the previous Sec- 
tion show that for low external magnetic fields (low vor- 
tex densities) the vortices move along approximately ID 
paths between the pillars. The effective ID potential 
landscape for the vortices has an approximately sinu- 
soidal form [see Eq. Q]. For simplicity we will consider 
here a periodic triangular potential such that the force 
on a vortex is given by: 

F(y) = -F P sign[cos(27ry/a)], 

which has the main ingredients needed to describe the 
soliton physics. Namely, spacial regions of weak Fc — Fp 
total force alternated with regions of large total force 
F c + F c . 

We also use a simpler form of vortex- vortex interaction 
potential with a cut-off: 



Here u = |r 12 /A, and the cut-off is at u c = 1. The 
precise form of the interaction is not important provided 



(v) sv = F c - 
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(Fc - Fp){F c + F P ) 
Fc 



(12) 



This result is valid for Fc > Fp, while for Fc < Fp the 
average velocity is zero (at zero temperature) and the 
dissipation vanishes. 

The increase in the slope of the magnetoresistance at 
the first matching field indicates an active role of the 
vortex- vortex interactions. To understand the behavior 
of the magnetoresistance at fields above the first match- 
ing field it is useful to calculate the velocity probability 
distribution for an isolated vortex 
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2F C 
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For Ap = Fc — Fp <C Fc there is a large probability of 
finding a vortex with a slow velocity Ap-P^In the process 
of going over a period of the potential, the vortex spends 
most of the time traveling at a slow velocity and it is the 
slow velocity which dominates the average velocity of the 
vortex. 

At the first matching field the vortices are evenly 
spaced and synchronized and their average velocity is 
simply given by Eq. (12). Adding a single vortex to 
the first matching field generates a soliton, i.e. a discom- 
mensuration that travels faster than the average velocity 
of the vortices. The extra vortex pushes its first neighbor 
which is in a slow velocity region increasing its velocity. 
The extra vortex stays trapped in a region with a slow 
velocity while its neighbor now with a fast velocity ap- 
proaches the next vortex in the chain to help it in turn 
go faster along the slow velocity region (see Fig. [7]). 

To illustrate this point we can consider a hard sphere 
potential of range a/2 for which it is easy to see that the 
extra vortex in the system amounts to the addition of a 
vortex with velocity v so i — Fc while the average velocity 
of the rest of the vortices is still given by the single vortex 
result. For the potential of Eq. ( 11 ) the average velocity 



of the soliton is slower than Fc but faster than the single 
vortex result. 

Increasing further the number of vortices increases the 
number of solitons and the dissipation. The increase is 
however larger than the one obtained for B < as the 
slope of R is given by the soliton velocity [see Eq. ^] . A 
weak soliton-soliton interaction leads to a linear increase 
of the resistivity. While the slope of the first segment 
is dominated by the single vortex result for the average 
velocity ~ Fc — Fp, the slope for fields above the first 
matching field is given by the speed of the soliton which 



FIG. 7. Schematic representation of a moving discommensu- 
ration in a ID system. The temporal sequence is starts at 
a) and ends at d). Initially there are as many vortices (filled 
red dots) as periods in the total potential (solid line) which is 
given by the pinning potential and the driving potential. A 
single vortex added to the system (open blue dot) generates 
a cascade of vortex displacements. 



is a fraction of Fc- For Fc ~ Fp, a small change in Fq 
can lead to a large relative change in the slope of the 
first segment while producing a small relative change in 
the slope of the second. This occurs in the experiments 
whenever there is a single kink in the magnetoresistance 
curve, which indicates that the velocity of the soliton is 
~ Fc and there is no room for a second increase in the 
slope at the second matching field. 

Increasing further the number of vortices up to the 
second matching field, assuming that the ID model still 
gives a qualitative picture for these vortex densities^! can 
lead to a second kink in the magnetoresistance (see Fig. 
[8]). The kink in the magnetoresistance at the second 
matching field is observed if the vortex-vortex interac- 
tion is weak enough (two vortices can be simultaneously 
in the same slow- velocity region). This can be observed 
in Fig. [8] where the two kink behavior with an associ- 
ated increase in the magnetoresistance is obtained. The 
dynamics above the second matching field is similar to 
the one above the first matching field. At the second 
matching field there are as many vortices in solitons as 
in the underlying vortex lattice. In order to observe a 
kink at 2B^, the soliton velocity must be <C Fq which 
means that the 'fast' vortex 1 pushing a 'slow' vortex 2 
remains trapped in the slow velocity region for some time 
before vortex 2 can leave it. In such situation, an extra 
vortex in the system can give an extra push to the pair 
of vortices accelerating the process and generating what 
we may dub a super-soliton with velocity V2 SO i- 

Decreasing the pinning force Fp or increasing the 
Lorentz force Fc reduces the time spent by the vortices 
in a slow velocity region and eventually destroys the pos- 
sibility of having a super-soliton regime. This is shown 
in the inset of Fig. [8] where we present simulations of the 
magnetoresistance for different values of Fc- 



0.5 1 1.5 2 2.5 

FIG. 8. Total velocity for vortices moving in a triangular 
potential. A two kink behavior is obtained as a function of the 
number of vortices. The dotted style lines are linear fits and 
serve as guides to the eye. The parameters are Fc = 1.03 Fp, 
a = 4, Fp — 10. The two kink behavior is lost for large enough 
values of Fc- This is shown in the inset for a different set 
of parameters F c = 1.025 F P , 1.075 F P , 1.275 -Fp, a = 4, 
Fp = 1.33. 



VII. SUMMARY AND CONCLUSIONS 

We have measured the electronic transport through pe- 
riodically modulated thickness Nb films in the supercon- 
ducting state as a function of the magnetic field. We 
observed a piecewise linear behavior of the resistance 
that we interpreted in terms of different dynamic regimes 
of the superconducting vortices. The thickness modula- 
tion of the film generates a periodic potential for the 
vortices which induces commensurability effects in the 
vortex flow. Numerical results for the vortex dynamics, 
as well as analytical results for a simplified model indi- 
cate the presence of different moving discommensuration 
(or soliton) regimes. In the regime of parameters where 
there is a good qualitative agreement with the experi- 
mental results, the solitons observed are well described 
considering independent one-dimensional vortex chains 
formed between pairs of neighbouring rows of cylinders. 
These results open the possibility of future studies that 
explore the crossover from a one-dimensional to a two- 
dimensional behavior of the vortex solitons. 

The generation of a thickness modulation using litho- 
graphic techniques offers a unique opportunity to con- 
trol precisely the intensity, geometry and topology of 
the pinning potential. Ginzburg-Landau calculations al- 
low in principle for a calculation of the pinning poten- 
tial surface, making these systems a fertile ground for 
experiment-theory comparisons and the study of differ- 
ent dynamical regimes and transitions. 
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Appendix A: ID sine potential 



For a cosenoidal one dimensional potential 



V(y) = a ^cos(2ny/a), 



the equation of motion of a single vortex 



(Al) 



y = F c - Fp sin(27ry/a) 



(A2) 



yit) 



tan 



F P 



+ (v) sv ( ir(v) sv t 7T 

— tan 

F c \ a 2 



(A3) 



F 2 
r p 



valid for —a/2 < y(i) < a/2, where (v) sv — \JF% 
is the average velocity of a vortex over a period of the 
pinning potential. The probability distribution of the 
vortex coordinate y over a period can be calculated using 
the relation P(y) — where P(t) — 



P(y) 



i 



a Fc — Fp sin(27ry/a) 



(A4) 



Using P(y) — ^fj- (here two branches must be added 
because the velocity is not a monotonous function of the 
coordinate) we finally obtain the single vortex velocity 
probability distribution 



Psv(V) 



J_ / (F c - F P )(Fc + F P ) 

irV V (V - F c + F P ){F C +F P -V) : 



(A5) 



which has square root divergences at U = Fc ± Fp . As in 
the toy model described in Sec. |VI[ for < Fc — F P <C 
Fc there is a large probability of finding the vortex with 
a low velocity V ~ Fc — Fp. This is an central ingredient 
for the soliton physics described in the main text. 
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